Paramagnetic Meissner effect in superconductors 
from self-consistent solution of Ginzburg— Landau equations 
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The paramagnetic Meissner effect (PME) is observed in small superconducting samples, and a num- 
ber of controversial explanations of this effect are proposed, but there is as yet no clear understanding 
of its nature. In the present paper PME is considered on the base of the Ginzburg-Landau theory 
(GL). The one-dimensional solutions are obtained in a model case of a long superconducting cylinder 
for different cylinder radii R, the GL-parameters X and vorticities m. Acording to GL-theory, PME 
is caused by the presence of vortices inside the sample. The superconducting current flows around 
the vortex to screeen the vortex internal field from the bulk of the sample. Another current flows 
at the boundary to screen the external field H from entering the sample. These screening currents 
flow in opposite directions and contribute with opposite signs to the total magnetic moment (or 
magnetization) of the sample. Depending on H , the total magnetization M may be either negative 
(diamagnetism) , or positive (paramagnetism). The detailed study of a very complicated saw-like 
dependence M(H) (and of other characteristics), which follow from the self-consistent solutions of 
the GL-equations, is presented. 



I. INTRODUCTION 

A superconductor, placed in magnetic field and cooled 
down through the transition temperature, expels the 
magnetic flux. This phenomen, known as the diamag- 
netic Meissner effect (DME), is one of the most essential 
properties of superconductors. Suprisingly, several re- 
cent experiments have shown that some superconducting 
samples may attract magnetic field, which corresponds to 
the so-called paramagneic Meissner effect (PME). Origi- 
nally, PME was observed in small size high-temperature 
superconductors (HTS), what has prompted the appear- 
ance of a number of theories, attributing this effect to 
a non-conventional superconductivity in these materials. 
(Mention such explanations, as spontaneous currents due 
to 7r-contacts, vortex pair fluctuations, d-wave supercon- 
ductivity, orbital glass, the Josephson junctions, the An- 
dreev reflection of pairs at the boundary, and others. The 
numerous references to the experimental and theoreti- 
cal work on PME may be found in recent papers [1-7].) 
Later, the anomalous PME was observed also in con- 
ventional low-T c supercontuctors and was considered as 
intrinsic property of any finite-size superconductor. 

In [2-5] the PME was studied in a case of conventional 
superconductors, using the Ginzburg-Landau (GL) sys- 
tem of equations. Here such conceptions were put for- 
ward, as the flux capture and its compression [2], the 
hysteresis transitions and the role of experimental setup 
[3] , the metastability of the vortex configuration due to 
the influence of the sample surface [4] , the disbalance of 
screening currents in metastable vortex states [5] . 

As was noticed in [5], the GL-theory treats PME in 
finite-size superconductors as purely electrodynamic ef- 



fect, caused by presence of vortices inside the supercon- 
ductor. If the vortex exists inside the sample, the super- 
conducting current j s flows around the vortex center in a 
direction to screen out the vortex field from entering the 
bulk of the sample. If the sample is placed in external 
magnetic field H , the additional surface current flows to 
screen out the field H from entering the sample. These 
two currents have opposite directions and produce con- 
tributions of opposite sign to the sample magnetic mo- 
ment (or, magnetization) M = (l/2c)/[j s x r]dv. The 
resulting value of M (in the direction of H z > 0) may be 
either negative (this corresponds to diamagnetic suscep- 
tibility), or positive (this corresponds to paramagnetic 
susceptibility), depending on H magnitude. 

In the present paper the origination of PME within the 
GL-thcory approach is studied in more details, as com- 
pared with [5], including both type-I and type- II super- 
conductors, for various sample dimensions and different 
values of GL-parameter x. In Sec. 2 the GL-equations 
for the order parameter ip and magnetic field B inside 
the superconductor are written. In Sec. 3 the results of 
numerical calculations of non-linear system of equations 
are presented by a number of graphics, accompanied by 
necessary comments. Sec. 4 contains a short resume and 
discussion of the results. 



II. EQUATIONS 

Below the case is considered of a long superconduct- 
ing cylinder of radius i?, in the external magnetic field 
H > 0, which is parallel to the cylinder element. In the 
cylindrical co-ordinates the system of GL-equations may 



be written in dimensionless form [5] 

tp 2 U = 0, 



d 2 U 
dp 2 



ldU 
P dp 



d 2 ip ldtp 2 
dp 2 p dp 



U 2 

^) - —^ = o. 
p 2 



(1) 

(2) 



Here U(p) is the dimensionless field potential; b(p) is the 
dimensionless magnetic field; ip(p) is the normalized or- 
der parameter; p — r/X, A is the field penetration length; 
A = x£, where £ is the coherence length, x is the GL- 
paramctcr. The dimensioned potential A, field B and 
current j s are related to the corresponding dimensionless 
quantities by the formulae [5]: 
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(The field B in (3) is normalized by H x = o /(2ttA 2 ), 
with b = B/H\: instead of H\ one can normalize by 
Hz = 0o/(2^ 2 ), or by H K = 0o/(2tt£A) = H s /x. The 
coefficients in (f), (2) would change accordingly.) The 
vorticity m in (3) specifics how many flux quanta are 
associated with the vortex, centered at the cylinder axis 
(the giant- vortex state). 

The boundary conditions to Eq. (1) are [5]: 
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where p\ = R/X, hx = H/H\. 

The boundary conditions to Eq. (2) are [5]: 
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The magnetic moment (or, magnetization) of the cylin- 
der, related to the unity volume, may be written in a form 



M 

V 



1 
v 



B-H 



4tt 



-dv 



B n 



H 



4n 



B av = 1 J B(r)dv = |*i, 

where B av is the mean field value inside the superconduc- 
tor, $i is the total magnetic flux, confined in the cylin- 
der. In the normalization (3), denoting b = B av /H\, 
h\ = H/H\, M\ = M/H\, one finds 
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Accordingly, the normalized Gibbs free energy of the 
system may be written as [5] 
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Here AG = G s — G n is the difference of free energies 
in superconducting and normal states; 6(0) = B(Q)/H\, 
B(0) is the magnetic field at the cylinder axis; H crn = 
0o/(2-7rv2A£) is thermodinamical critical field of massive 
superconductor; go is the condensation energy with ac- 
count for the order parameter gradient. The expressions 
(6), (7) will be used below for calculating the correspond- 
ing quantities. 



III. NUMERICAL RESULTS 

The system of equations (f), (2) was solved numeri- 
cally by the method, described in [5,8]. For every fixed 
set of parameters R, H, to, x the unique self-consistent 
solution for ip(p) and U(p), satisfying the boundary con- 
ditions (4), (5), was found. If the external field H exceeds 
some critical value He (H > He ), no superconduct- 
ing solution with ijj ^ is possible. The phase diagrams, 
which separate the superconducting (ip ^ 0) and normal 
(ijj = 0) regions, are presented in Fig. 1 for different 
vorticities to and parameters x. 

The critical field of the vortex-free Meissner state 
(to = 0) is distinguished by thick curve. The super- 
conducting region lays between this curve and the axis 
/if = (/i^ = H/H^). For to > the superconducting 
region is inside of the corresponding m-curvc (thin line) . 
To every point of the superconducting TO-region corre- 
sponds some solution of Eqs. (1)— (5). The transforma- 
tion of solutions, when the superconducting regions are 
crossed along the lines Rx = R/X = const (see the dotted 
horizontal lines), is illustrated below by Figs. (2)-(4). 

In Figs. l(a, 6) the dotted vertical line h^ — 1.695 
represents the critical field H C 3 = 1.695i/ C 2 (H C 2 = 
0/(27r£ 2 ) = H^), at which the surface superconductiv- 
ity [9,10] nucleates in large (Rx 3> 1) samples. It is seen, 
that in type-II superconductors (x > l/\/2) the super- 
conducting states (m > 0) are possible, which lay to the 
right of the line m = 0. One may say, that these states 
are enhanced by the external field. In contrast to this, in 
type-I superconductors (see Fig. 1(c)) the states to > 



lay to the left of the line to = 0, so they are weakend by 
the field. 

Fig. 2 illustrates the superconductor behavior ( R\ = 
2, x = 2), when the field h^ varies. In Fig. 2(a) the re- 
gions are depicted, where the state m can exist (ip m ax is 
the maximal value of the order parameter in the corre- 
sponding m-state). For instance, the state to — 4 (which 
is distinguished for clarity by thick line) can exist only in 
the interval of fields, where ip m ax > 0. The state to = 11 
can exist only in sufficiently large fields, h% ~ 2 (see also 
Fig. 1(a)). 

Fig. 2(6) shows the field dependence of the cylinder 
magnetization in different m-states. For instance, the 
magnetization (— 4ttM\) in the state to = 4 may be either 
positive (this corresponds to diamagnetic susceptibility, 
M\ < 0), or negative (this corresponds to paramagnetic 
susceptibility, M\ > 0), depending on the field value. 
The vortex-free Meissner state (to = 0) is diamagnetic 
for all fields. 

Fig. 2(c) shows the normalized free-energy difference, 
Ag, for the states with different to. At the points, where 
the curves with neighboring values of to cross each other 
(marked by open circles), the equilibrium transitions be- 
tween the states to and to ± 1 became possible. These 
equilibrium transitions are marked by the dashed vertical 
lines in Fig. 2(6). If the equilibrium transitions are only 
allowed, the magnetization curve would always be dia- 
magnetic (— A-kM\ > 0), and have a reversible saw- like 
shape. However, if the metastable states are also allowed, 
the magnetization curve would look very different. This 
is illustrated by Figs. 3(a, 6, c) for x = 1 and by Figs. 
3(d, e, /) for x = 0.5. 

Consider Fig. 3(c), the state m = 2. That part of 
the curve to = 2, which lays to the left of the equilibrium 
point (eg), is metastable (because the lower energy states 
to = 0, or m — 1, exist). If the field h^ is diminished from 
its equilibrium value (in the field-cooling regime), the 
state to = 2 can persist only down to the point a, where 
the forced transition to other state should occure. The 
nearest laying energy state is to = 1, it corresponds to 
the relative minimum of the free energy (the state m = 
corresponds to the absolute minimum) . If the transitions 
between the nearest laying energy minimae are possible, 
then at the point a the transition from the state m — 2 
to the state m = 1 may occure. In Fg. 3(6) this would 
mean the transition from the paramagnetic state to = 2 
(which lays to the left of the dashed-dotted vertical line, 
passing through the point —4irM\ = 0) to the paramag- 
netic state to = 1. If the field h^ is now increased, the 
state to = 1 may survive up to the point S, where the 
forced transition from the diamagnetic state to = 1 to 
the diamagnetic state m = 3 may occure. At the point (5 
the transition from paramagnetic part of the curve m = 3 
to the paramagnetic part of the curve to = 2 occures. At 
the point 7 the transition between diamagnetic (to = 0) 
and paramagnetic (to = 3) states may happen. Note, 



that all paramagnetic states are metastable (they lay to 
the left of the corresponding equilibrium points). Evi- 
dently, the transitions between various states, shown in 
Figs. 3(a,6, c), may lead to the hysteresis behavior of 
M(H) (see also [3,4]). 

[Note, that for larger values of R\ and x the reversible 
jump transitions to the edge-suppressed states [5] of the 
same vorticity to are also possible, which might addition- 
ally complicate the dependence M(H).] 

The analogous transitions between paramagnetic and 
diamagnetic states are also possible in type-I supercon- 
ductors (see Figs. 3(d, e, /)). Notice, that the curve 
to = terminates at the point 7, where Ag > 0. This 
part of the curve corresponds to the "superheated" su- 
perconducting state, considered first by Ginzburg [11] for 
type-I superconductors in the approximation ip = const 
(see also [12]). The point 7 corresponds to the maximal 
superheating. The only possible transition here may be 
to the normal state (ip — 0). If the field h^ is now di- 
minished, the normal state may be "supercooled" [11,12] 
down to the point 5 (fi£ — 1.49), where the supercon- 
ducting state recovers. The arrows in Fig. 3(d) mark the 
boundary of the hysteresis loop. (The supercooled states 
exist also in a case x > l/\/2- They additionally com- 
plicate the picture and are not shown in Figs. 3(a, 6, c).) 

As is seen from Figs. 2,3, each magnetization curve 
(— 47rM^ ) has the diamagnetic and paramagnetic parts. 
Consider points of maximal diamagnetism and maximal 
paramagnetism. Fig. 4 illustrates the space distribu- 
tions i[>(r), 6(r), j(r) at these points for several val- 
ues of m,x,R\. The curves 1 (thick lines in Fig. 4) 
correspond to maximal paramagnetism (M\ > 0), the 
curves 3 (thin lines) correspond to maximal diamag- 
netism (M\ < 0), the curves 2 (dashed lines) represent 
points, where M\ = 0. 

Figs. 4(6, e, h) show, that for paramagnetic curves 1 
the magnetic field inside the sample is greater, than the 
applied field h\ (6(r) > h\, h\ = x 2 h^). This means, 
that the influence of the internal vortex field -B(O) is 
greater, than the influence of the external field H . The 
relative suppression of the order parameter, seen for the 
curves 1 in Figs. 4(a, d, g), is caused mainly by the inter- 
nal field -B(O). The same is seen in Figs. 4(c, /, i), where 
all the currents 1 are paramagnetic (j > 0), they screen 
mainly the internal vortex field. 

For larger fields H (the curves 2) , the role of external 
magnetic field is more pronounced. Now, partly 6(r) > 
h\ (the paramagnetic behavior), and partly b(r) < h\ 
(the diamagnetic behavior). Analogously behaves the 
current j: part of the current is j > (paramagnetic 
screening of the vortex field), and partly j < (diamag- 
netic screening of the external field). These two parts of 
j counterbalance each other, resulting in zero magneti- 
zation, M\ = 0. 

The order parameter tp(r) in the case 2 is enhanced by 



a larger field (compare the curves 2 and 1 for ip). This 
reflects the reentrant behavior of small-size samples: the 
superconductivity, which was impossible in small exter- 
nal field, is possible again in larger fields (see Figs. 2,3). 

When the field H is increased further (the curves 
3 in Fig. 4), the sample passes to diamagnctic state 
(M\ < 0): the field b(r) < h\ everywhere, the current j 
screens predominantly the external field (j(r) < near 
the surface), the order parameter ip(r) displays the sup- 
pression by strong external field H . 

In the foregoing discussion the case was considered, 
when the external field H has been varied, but the tem- 
perature T was fixed (see the horizontal dotted lines 
T = const in Fig. 1). In some experiments (for in- 
stance, [7]) the field H was fixed, but the temperature 
T varied. In this case the transitions between different 
diamagnetic and paramagnetic states (see Figs. 3(6, c)) 
are also possible. 

Indeed, in the vicinity of T c the parameters £ and 
A = x£ depend on the sample temperature T by the law 
£(T) = £ (1 - T/Tc)- 1 / 2 . Using notations: t = T/T c , 
R x = R/X, h£ = H/Hfr H^ = </> /(2tt£ 2 ), one has two 
relations: R\ = ro\/I — t and h^ = ho/(l — t), or 
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where r = R/X , h = H/H i0 , H i0 = </> /(27r£ 2 ), A = 
x£o- From Eqs. (8) the dependence between R\ and h^ 
follows: 
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The dashed curves in Fig. 1(6) are the lines, along which 
the representation point (R\, h^) moves, when H = const 
and the temperature T varies only. Taking x = 1, 
H = 100 Oe, one has A = 1 for R = 1.8 ■ 10 _5 cm, 
and A = 3 for R = 5.4 • 10~ 5 cm. Notice the essential 
difference in the dashed curves positions. For smaller 
R (A = 1) the system remains mainly in the diamag- 
netic state to — 0, so no PME is possible. For larger R 
(A = 3) the dashed curve H = const crosses several of re- 
state phase boundaries. This more complicated situation 
is illustrated by Fig. 5. 

In Fig. 5 the temperature dependencies are shown: 



of the order parameter ip 



(m) 



of the magnetization 



(— AttM^ 1 '), and the Gibbs free energy A</ m ) in several 
m-states. It is seen, that the first superconducting state 
(ip ^ 0), encountered in the temperature-cooled mode, is 
m = 2 (see Fig. 1(6), the point /i). This is the state of 
the surface superconductivity in a finite-dimension sam- 
pie. The magnetization curve is diamagnetic (M^ < 0), 
it starts from zero, as the second order phase transition. 
When the temperature t — T/T c is reduced further, other 
m-states became possible. If the temperature t is swept 



up and down (as in [7]), alongside with the equilibrium 
transitions between diamagnetic states (at the points a 
and /? in Figs. 5(6, c)) the transitions between various 
metastable branches of the curves may also occure. For 
instance, at the termination point 7 the jump transi- 
tion from the diamagnetic state m = to the nearest 
energy laying paramagnetic state m — 4 is possible; at 
the termination point 8 the jump from the paramagnetic 
state to = 4 to the diamagnetic state m = 1 may oc- 
cure. The up and down arrows at the points 7 and S in 
Fig. 5(6) mark the possible hysteresis loop in the M(T)- 
dependence (compare with the experimentally found de- 
pendencies [7]). 

A more detailed comparison of the present calculations 
with experiments on PME is not possible, because the 
model case of a long circular cylinder approximates the 
real specimens geometry very crudely. However, it is in- 
structive to see the results of GL-theory, because they 
suggest some ideas, which should be viewed with atten- 
tion. For instance, in the scope of GL-theory the PME 
looks more like a bulk effect, it is caused by the balance of 
currents, flowing inside the sample. The role of the sur- 
face seems not to be crucial, though in the finitc-hight 
extended samples [7] the role of the boundary may be 
important. It also looks probable, that PME might be 
observed in the proximity structures [6] , where the super- 
conducting cylinders are embodied into the normal-metal 
enviroment, or in a system of microholes inside the su- 
perconducting matrix [13]. 

Below are comments on some of the papers, where the 
PME was first considered within the GL-approach. 

In [15] the hypothesis was put forward of a normal 
layer, which might be formed inside the superconducting 
plate under the influence of the external field H and is 
situated at some distance from the plate boundary. Ac- 
cording to [15], the conception of such a layer might help 
in understanding the nature of PME in superconductors. 
Note, that in the homogeneous GL-theory no such layer 
is possible, though the basic idea of [15] is, probably, 
correct. Indeed, if a chain of vortices is formed near the 
plate surface (with the order parameter vanishing at their 
axes), the average vortices distribution may be visualized 
as a normal layer, in accordance with [15]. 

In [2] the conception of the flux compression was pro- 
posed to explain the PME. According to [2], the giant 
TO-vortex state, which is formed at sufficiently high field 
H, persist at smaller H in a form of the flux-compressed 
state, with the eventual onset of the paramagnetic re- 
sponse. Some evidence of the possible PME in a small- 
radius cylinder was obtained in [2] by self-consistent so- 
lution of GL-cquations. Note, that our calculations are 
along the same lines (see Figs. 2, 3 and also [5]), but 
they show, that the giant m-vortex state can persist only 
in a finite field interval, outside which the vorticity to 
changes. As it follows from Figs. 4(6, e, ft-), the impo- 
sition of the external field H leads not to a field com- 



pression, but, better to say, to a field submission effect. 
Indeed, at smaller H (the curves 1) the external field is 
effectively screened out (b(r) > h\); at larger H (the 
curves 2) the Meissner screening is not so effective and 
the external field is partly submitted into the sample 
(h\ > b(r) near the boundary); still at larger H (the 
curves 3) the external field floods almost all the interior 
(h\ > b(r) everywhere), with eventually total submission 
to the external field (b(r) = h\ in the normal state). As 
was mentioned above, the appearance of PME in this pic- 
ture is due to the redistribution of the screening currents, 
flowing inside the specimen. 

[There exist some doubts concerning the self- 
consistency of the solutions, found numerically in [2]. 
Thus, in Fig. 4(a) [2] the order parameter must depend 
linearly on r for r <C 1. The local induction b(r) in 
Figs. 4(c) and 5(a) [2] is maximal not at r = (where 
the vortex axis is situated), but at some distance r =/= 0, 
which seems strange, etc. Most closely our calculations 
correlate with those in [3,4,15], where many fragments 
of the picture, presented above, may be found, though 
no concise formulation of a possible nature of PME in 
conventional superconductors was given.] 



IV. CONCLUSIONS AND DISCUSSION 

The GL-theory of a long superconducting cylinder, 
which accounts for the transitions between centrally sym- 
metric m-states (the giant- vortex states), gives only qual- 
itative picture of the processes inside the real supercon- 
ductor. In reality, the giant-vortex may split itself into 
the asymmetric multy- vortex configuration of the same 
vorticity m. The specimen finite hight in z-direction, 
and the boundary conditions at the surface, might also 
be important for calculating the real vortex structure. 
All this makes the analysis more difficult. [The approx- 
imate GL-description of the multy-vortex configurations 
in mcsoscopic 3-dimensional samples may be found in 
[3,4] .) However, even in a general case of arbitrary vortex 
configuration, there will be two currents, flowing in oppo- 
site directions: one - to screen out the combined fields of 
the vortices, and second - to screen out the external field. 
Thus, basing on the GL-model calculations, it seems rea- 
sonable to offer the following qualitative scenario for the 
appearance of PME in finite-dimension samples. 

According to GL-thcory, the paramagnetism in con- 
ventional low-Tc superconductors may be caused by the 
presence of vortices inside the superconductor. The vor- 
tices, which are formed inside the sample in the field- 
cooling regime, may be pinned to the boundary, which 
is a source of inhomogeneity The directions of the ex- 
ternal and of the vortex magnetic fields are the same, 
but the directions of currents, which flow to screen these 
fields from the bulk of the sample, are opposite (clock- 
wise and counter clock-wise). The screening currents are 



responsible for two contributions of opposite sign to the 
total magnetic moment (or, equivalently, magnetization) 
of the sample: M — M para + Mdia- The distributions 
of the fields and currents inside the sample varies with 
H , so the measured magnetization may be either positive 
(M > 0, PME), or negative (M < 0, DME), or M = 0, 
depending on H. 

The pinning of the vortices to the boundary is a non- 
linear effect, described by the system of GL-equations. 
Each of the vortex configurations of fixed vorticity m can 
exist only within a limited interval of fields H . When the 
field takes the value outside this interval, the existing vor- 
tex configuration reconstructs, and a new configuration 
with different vorticity to' establishes (see Fig. 3). The 
jump on the magnetization curve may accompany such 
reconstruction. The new configuration may reveal cither 
PME, or DME, depending on the resulting state char- 
acteristics. If the system may occupy only the ground 
state of minimal free energy, then no PME is possible, 
and the behavior M(H) should be completely reversible. 
The presence of PME and of the hysteresis is a signal of 
metastability The examples of transitions between var- 
ious metastable states are presented in Fig. 3, both for 
type-I and type-II superconductors. 

The proposed qualitative GL-scenario looks rather 
general and one may attempt to apply it also to the case 
of HTS (x» 1), where PME is also observed in a field- 
cooling regime, is also characterized by the hysteresis, by 
the jump transitions between PME and DME branches 
of magnetization, and by the transitions between various 
quantum flux states (see, for instance, [14]). Whether 
PME in HTS can be explained by the presence of vor- 
tices, or by more exotic mechanisms, mentioned in the 
Introduction, is a question for further experimental and 
theoretical study. 
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x 2 ^): l-h\ = 1.52, 2-h\ = 2.921, 3 - h\ = 4.30 (if 
R x = 2, x = 2, m = 4); 1 - h\ = 0.45, 2 - hx = 0.985, 3 

hx = 1-52 (if R x = 2.5, x = 1, m = 2); 1 - hx = 0, 2 - 
h x = 0.195, 3-h x = 0.48 (if R x = 4, x = 0.5, m = 1). 

Fig. 5. Shown are as functions of temperature t = 
T/T c : (a) the maximal value of the order parameter, 
ipmax', (b) - the magnetization (— 4ttMx); (c) - the nor- 
malized free-energy difference, Ag. [For m = 0,1,2,3,4 



and A = 3 (x = 1, R 

Fig. 1(6) and Eq. (9).] 



5.4- 10~ 5 cm, H = 100 Oe), see 



Figures captions 

Fig. 1. The superconducting phase boundaries of vor- 
tex states with different vorticities m (the numerals at 
the curves), (a) - x = 2; (6) - x = 1; (c) - x = 0.5. The 
vertical dotted line h% = 1 corresponds to the critical field 
H C 2, the line h% — 1.695 corresponds to the critical field 
H c3 . (Here R x = R/X, ^ = H/H ( , H^ = 0o/(2< 2 ).) 

Fig. 2. (a) - The maximal value of the order parame- 
ter, ipmax, in m-state; (b) - the magnetization (— 47rM.\), 
Eq. (6); Mx = M/H x , H x - ^o/(2tt\ 2 ); (c) - the 
normalized free-energy, Eq. (7), as functions of h% for 
Rx = 2, x=2. 

Fig. 3. The same, as in Fig. 2, but for: (a, b, c) - 
Rx = 2.5, x = 1; (d, e,f)-R x = 4,x= 0.5. The letters 
a, /?, 7, (5, e mark the points, where m-solution terminates; 
the transitions to the nearest laying energy states are 
shown by the vertical arrows. The paramagnetic parts of 
the curves lay to the left of vertical lines, where Mx — 
0. [The arrows in Fig. 3(a) show the transition to the 
normal state (-0 = 0). The arrows in Figs. 3(d) mark the 
possible hysteresis loop boundaries for the state m = 0: 
7 - superheated (sh) Meissner state, 5 - supercooled (sc) 
normal state.] 

Fig. 4. The co-ordinate dependences: (a, d, g) - the or- 
der parameter -0, (6, e, h) - the magnetic field 6, (c, /, i) 
- the current j; for solutions with different Rx,x,m (see 
numbers in the insets) . The curves 1 correspond to max- 
imal paramagnetism (M > 0), the curves 2 correspond 
to M = 0, the curves 3 correspond to maximal diamag- 
netism (M < 0). The curves 1,2,3 are calculated for the 
following values of the applied field hx = H/Hx (hx = 
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